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Abstract Observational data require a rich variety of mechanisms to accelerate fast 
particles in astrophysical environments operating under different conditions. The 
mechanisms discussed in the literature include varying magnetic fields in compact 
sources, stochastic processes in turbulent environments, and acceleration behind 
shocks. An alternative, much less explored mechanism involves particle accelera- 
tion within magnetic reconnection sites. In this chapter we discuss this mechanism 
and show that particles can be efficiently accelerated by reconnection through a first 
order Fermi process within large scale current sheets (specially when in the presence 
of local turbulence which speeds up the reconnection and make the acceleration re- 
gion thicker) and also through a second order Fermi process in pure MHD turbulent 
environments. 



1 Introduction 

Energetic particles are ubiquitous in astrophysical environments and their accelera- 
tion still challenges the researchers. For instance, the origin of the ultra high energy 
cosmic rays (UHECRs) is unknown yet. Their spectrum is consistent with an origin 
in extragalactic astrophysical sources and candidates range from the birth of com- 
pact objects to explosions related to gamma-ray bursts (GRBs), or events in active 
galaxies (AGNs) ll2Tll . however, the mechanism(s) that produce(s) them is(are) still 
not fully understood. Similarly, recent very high energy observations with the Fermi 
and Swift satellites and ground based gamma ray observatories (HESS, VERITAS 
and MAGIC) of AGNs and GRBs have been challenging the current particle accel- 
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eration theories which have to explain how particles are accelerated to TeV or larger 
energies in regions relatively small compared to the fiducial scale of their sources 

ma. 

The mechanisms frequently discussed in the literature for accelerating ener- 
getic particles include varying magnetic fields in compact sources (e.g., Il2l [31 1411 '). 
stochastic processes in turbulent environments [41], and acceleration behind shocks. 
The latter, in particular, has been extensively discussed in the literature [48, 411 E] 
(see also the Chapter by Jokipii & Giacolone in this book). An alternative, much 
less explored mechanism so far, involves particle acceleration within magnetic re- 
connection sites. 

Magnetic reconnection occurs when two magnetic fluxes of opposite polarity 
encounter each other In the presence of finite magnetic resistivity, the converging 
magnetic field lines annihilate at the discontinuity surface and a current sheet forms 
there. 

In 2005, de Gouveia Dal Pino & Lazarian [4 [(henceforth GL05) proposed a 
mechanism to accelerate particles to relativistic velocities within the reconnection 
zone in a similar way to the first-order Fermi process that occurs in shocks. It is 
known from shock acceleration theory that particles are injected upstream and al- 
lowed to convect into the shock, while diffusing in space so as to effect multiple 
shock crossings, and thereby gain energy through the first order Fermi process [ 16 1. 
Similarly, GL05 [4| proposed that trapped charged particles may bounce back and 
forth several times and gain energy due to head-on collisions with the two converg- 
ing magnetic fluxes of opposite polarity that move to each other at the reconnec- 
tion velocity (Vr). They found that the particle energy gain after each round trip is 
AE/E o^Vr/c. Under fast magnetic reconnection conditions, e.g. induced by turbu- 
lence ll28l . Vr can be of the order of the local Alfven speed Va (see below and also 
the Chapter by Lazarian in this book. At the surroundings of relativistic sources, 
for instance, Vrec ~ ~ c, so that the mechanism can be rather efficient. GL05 [4 | 
have also shown that the accelerated particles have a power-law distribution and a 
corresponding electron synchrotron radio power-law spectrum which is compatible 
with the observed radio flares of galactic microquasars. 

Afterwards, Drake et al. (2006) [10] invoked a similar process, but within a col- 
lisionless reconnection scenario. In their model, the contraction of two-dimensional 
magnetic loops is controlled by the firehose instability that arises in a particle-in-cell 
(PIC) domain 140) IT2l . Other processes of acceleration, e.g. due to the electric field 
associated with the current in the reconnection region [37j and turbulence arising as 
a result of reconnection [27] were shown to be less dominant. 

Magnetic reconnection is very frequent and therefore it should be expected to 
induce acceleration of particles in a wide range of galactic and extragalactic envi- 
ronments. Originally discussed predominantly in the context of electrons in solar 
flares ilTOl im ITSl [19] [541 . it was later applied to explain the origin of anomalous 
cosmic ray protons [[311 [T2l . and the anisotropics in the direction of solar system 
magnetotail |[32]| . It also has been gaining importance in more extreme astrophysi- 
cal environments and sources, such as in the production of ultra high energy cosmic 
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rays ||2TI . in particle acceleration in jet-accretion disk systems |l4l|5]|6][I3[8l, and in 
the general framework of compact sources, as and GRBs ||29] |6] [TTl l53l 151] 1521 171 . 

The applications above, however, still require deep studies of the particle accel- 
eration in magnetic reconnection sites, as well as its connection with magnetohy- 
drodynamical (MHD) turbulence and fast magnetic reconnection. 

Steps in this direction has been performed recently fS] H ED |23] |24l where the 
model of GL05 [4] was tested successfully by means of two (2D) and three di- 
mensional (3D) MHD simulations. Kowal, de Gouveia Dal Pino & Lazarian (2011, 
henceforth LGKll) |23|, in particular, have shown that the acceleration of particles 
inserted in MHD domains of reconnection without including kinetic effects pro- 
duces results similar to those found in particle-in-cell (PIC) (collisionless) simula- 
tions where particle acceleration is controlled by kinetic effects such as the firehose 
instability |fT0l[T2l . This demonstrated that the acceleration in reconnection regions 
is a universal process which is not determined by the details of the plasma physics 
and can be also very efficient in coUisional gas. They have also shown that particle 
acceleration in 3D MHD reconnection behaves quite differently from the acceler- 
ation in 2D domains because the increase in the acceleration component parallel 
to the magnetic field is not constrained by the production or size of contracting is- 
lands, as in the 2D case. These results call for focusing on realistic 3D geometries 
of reconnection. 

Other concomitant studies have also explored test particle acceleration in MHD 
domains 1 18, 19 1. Gordovskyy et al. (2010) 1 18 |, for instance, focused on 2D models 
with time-dependent reconnection, while Gordovskyy and Browning (2011) llT9l . 
aiming the study of particle acceleration in small solar flares, examined a somewhat 
different scenario, with the acceleration of test particles by magnetic reconnection 
induced by kink instabilities in 3D twisted magnetic loops. Although they have ob- 
tained results for the particle energy distributions which are compatible with field- 
aligned acceleration as in the studies above, they did not explore the nature of the 
mechanism accelerating the particles. 

Kowal, de Gouveia Dal Pino & Lazarian (2012, henceforth KGL12)|24| have, 
in turn, injected test particles in different 3D MHD of reconnection and compared 
the particle spectrum and acceleration rates in these different domains. When con- 
sidering a single Sweet-Parker topology ll50l l44l (subject to large artificial magnetic 
resistivity to allow fast reconnection), they have found that particles accelerate pre- 
dominantly through a first-order Fermi process, as predicted in GL05 [4J. When 
turbulence is included within the current sheet, the acceleration is highly enhanced. 
This is due to the fact that reconnection becomes fast in a natural way (and inde- 
pendent of magnetic resistivity) in the presence of turbulence Il28l l22l and allows 
the formation of a thick volume filled with multiple simultaneously reconnecting 
magnetic fluxes. The particles trapped within this volume suffer several head-on 
scatterings with the contracting magnetic fluctuations, which significantly increase 
the acceleration rate and the amount of particles which are accelerated through a 
first-order Fermi process. They have also tested the acceleration of particles in a 
pure MHD turbulence, where particles suffer collisions both with approaching and 
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Fig. 1 Particle acceleration in a reconnection site where two magnetic field fluxes of opposite 
polarity move to each other. Cosmic rays spiral about a reconnected magnetic field line and bounce 
back and forth at points A and B within the current sheet. The reconnected regions move towards 
each other with the reconnection velocity Vg. Particles gain energy due to "collisions" with the 
magnetic iiregularities within the two converging fluxes, j ust like in the first-order Fermi process 
in shock fronts (see GL05 |4|). The advection of cosmic rays entrained on magnetic field lines 
happens at the outflow velocity, which in most cases is of the order of Va- Bouncing between 
the points A and B happens either because of streaming instability induced by energetic particles 
or by magnetic turbulence in the reconnection region (as discussed in |28|). In particular, when 
turbulence is present, the acceleration region is filled in by several oppositely moving reconnected 
flux tubes which collide and repeat on smaller and smaller scales the pattern of the larger scale 
reconnection making the process very fast and therefore, the particle acceleration very efficient. 
From 1331 . 



receding magnetic irregularities. The acceleration rate is smaller in this case and 
suggests that the dominant process is a second order Fermi. 

In this chapter, we discuss these acceleration mechanisms in magnetic recon- 
nection sites in detail and review the recent analytical and numerical results in this 
regard. 



2 An analytical model for first order Fermi particle acceleration 
within magnetic reconnection sites 

We first discuss an analytical model for acceleration of particles within reconnection 
sites which was originally introduced by GL05 (U (see also ||30l l5l. Figure [T]illus- 
trates the simplest realization of the acceleration within a large scale reconnection 
region. As described in the Figure, as a particle bounces back and forth between two 
converging magnetic fluxes of opposite polarity, it gains energy through a first-order 
Fermi acceleration. 

In order to derive the energy spectrum of the accelerated particles one can in- 
voke a similar procedure to the one employed in the calculation of the first order 
Fermi acceleration in shocks (see, e.g., |[T] |39] |33]). Let us consider the acceler- 
ation of Mo particles with an initial energy Eq. If a particle acquires an energy 
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E = PEq after a collision, its energy after m collisions will I5'"Eq. At the same 
time if the probability of a particle to remain within the acceleration region is P, 
after m collisions the number of accelerated particles will be M = P"'Mo. Therefore, 
ln(M/Mo)/ln(£'/£'o) =lnP/lnj3 or 

Mo " Uo j 

Since some of these M particles will be further accelerated before escaping the sys- 
tem, the equation above implies that the number A^(£') of particles accelerated to 
energies equal to or larger than E is given by: 

N{E)dE = const x (2) 

To compute P and j3 within the reconnection site we may consider the follow- 
ing process. The particles from the upper reconnection region will "see" the lower 
reconnection region moving towards them with the velocity 2Vr (see Figure [U. 
If a particle from the upper region enters at an angle 9 with respect to the direc- 
tion of Vr into the lower region then the expected energy gain of the particle is 
5E/E = 2Vr cos jc. For an isotropic distribution of particles their probability func- 
tion is = 2 sin cos QdQ and therefore the average energy gain per crossing of 
the reconnection region is 

Vr I'^I^ n Vr 

{8E/E) = — 2cos-0sin0t/0 =4/3— (3) 
c Jo c 

Particles will complete a full acceleration cycle when they return back to the upper 
reconnection region. Similarly, if they are in the lower reconnection region they 
will see the upper reconnection region moving towards them with the speed 2Vr. 
As a result, a full acceleration cycle provides an energy increase {5E /E)cYcie = 
8/3 (Vr/c) and thus 

I5^E/Eo^I+8/3{Vr/c) (4) 

Let us assume that the particle diffusion velocity is much smaller than Vr. In 
analogy to particle acceleration in a shock front, for simplicity, we further assume 
that the total number of particles crossing the boundaries of the upper and lower 
magnetic fluxes is 2 x I /4(nc), where n is the number density of particles. If the 
particles are advected to outside of the reconnection region with the magnetized 
plasma outflow then, the loss of the energetic particles will be given approximately 
by IVrh. Therefore the fraction of energetic particles which are lost in a cycle will 
be y«n/[l /4(«c)] = 4Vr/c and 

P^I-Wr/c. (5) 
Combining Eqs. (|4]i, and ^ one obtains 



N{E)dE = constiE-^l^dE, 



(6) 
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which is the spectrum of accelerated energetic particles for the case when the plasma 
back-reaction is negligible (|4]|33]. We note that the power-law index obtained above 
is independent of the reconnection velocity Vr. This is in part due to the simplified 
assumptions in the derivation above. Nonetheless, we will see below that in the 
much more realistic numerical simulations of the acceleration of test particles in 
non-relativistic MHD reconnection sites, the acceleration rate and power spectrum 
are not much sensitive to the reconnection speed [9 |. 

In recent work, Drury (2012) lfl4l tried to improve the analytical model above by 
considering two additional effects. First, he took into account the energy losses due 
to the outflow from the reconnection region (which is associated with a divergence 
of the flow field) and second, he relaxed the assumption considered above (GL05) 
that the escape rate is the same as that from a shock. Then, he repeated the calcula- 
tion above and obtained a power law spectral index which is the same as in shock 
acceleration if expressed in terms of the compression ratio in the system r — P2/P1, 
where pi and p2 are the plasma densities at the inflow and the outflow regions of 
the reconnection site, respectively. In other words, he obtained; 

f{p)^P^, (7) 

where f{p) is the particles distribution function as a function of their momentum p. 
Considering that only little energy is used to heat the plasma and the conversion is 
essentially one of magnetic energy into kinetic energy, then the reconnection can be 
very compressive, unless the outflow is significantly over pressured relative to the 
environment. But, for a strongly magnetized inflow this is a very weak constraint 
lfT4ll . Thus one can expect in general large values of the compression ratio, possibly 
larger than the value four usually assumed for adiabatic shocks. For large values of 
r one finds f{p) °^ p^^, or: 

N{E)ocE-^ (8) 

which is a power-law spectrum much harder than the one obtained above by GL05, 
but confirms their prediction that a rather efficient first-order Fermi particle acceler- 
ation process can take place in magnetic reconnection sites. 

The considerations above also allow one to estimate the acceleration time-scale 
due to reconnection, which is a straightforward generalization of the result for shock 
acceleration as well. If the accelerated particles penetrate a distance Li into the 
inflow and L2 into the outflow region, then one can show that the acceleration time 
scale is approximately given by lfT4l 

A1L1+A2L2 + W 
^^"■^ A1U1-A2U2 ^' 

where W is the volume of the reconnection region; Ui and A i are the velocity and 
the cross-sectional area of the reconnection region for the inflow, respectively; and 
U2 and A2 are the same for the outflow. This acceleration time scale is similar to 
that for shock acceleration, and the constraints on maximum energy due to the finite 
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age and size of the reconnection region will thus be comparable to those in shock 
acceleration. 

It should be noticed that Drury's model [ 14] predicts that the acceleration within 
reconnection sites requires a large compression ratio in order to be efficient. How- 
ever, according to the discussion in the previous section (see also the following 
sections), the requirement for the magnetic reconnection acceleration process to 
proceed efficiently is to keep the accelerated particles within contracting magnetic 
loops. This requires constraints on the particle diffusivity specially perpendicular to 
the magnetic field direction. The subtlety of this point is related to the fact that while 
in the first-order Fermi acceleration in shocks magnetic compression is important, 
the acceleration via reconnection is applicable also to incompressible fluids. Thus, 
unlike shocks, not the entire volume has to shrink for the acceleration to occur, 
but only the volume of the magnetic flux tube. If the perpendicular diffusion of the 
particles to the magnetic field is large they may decouple from the magnetic field. 
Indeed, it is easy to see that as long as the particles in the magnetic flux rope, as de- 
picted in Figure 1, bounce back and forth between the converging mirrors they are 
accelerated. However, if these particles leave the flux rope too fast, they may start 
bouncing between the magnetic fields of different flux ropes which may sometimes 
decrease their energy. Thus it is important that the particle diffusion parallel and 
perpendicular to the magnetic field stays different. Particle anisotropy which arises 
from particle preferentially being accelerated in the parallel direction may also be 
important. 

In the case of laminar reconnection sites (i.e., with no turbulence), as in a Sweet- 
Parker model for reconnection (and as assumed above; also see Section 3.4), the 
scales can be much smaller compared to shock structures. This will lead to lower 
maximum energies (the dominant limiting loss process being diffusion out of the 
sides of the reconnection region). This means that the models developed above ap- 
ply only if the particle diffusion length scale in the inflow is small compared to its 
lateral extent. Thus they may only apply in a restricted energy range which should 
be contrasted with the case of shock acceleration where the scale separation is, in 
general, much larger However, when turbulence is present within the reconnection 
site, this will make the reconnection volume much larger Il28ll22ll23, 24 1 and there- 
fore, this process can be competitive to shock acceleration, as we will see below. 



3 Particle Acceleration in reconnection sites: numerical studies 

As discussed in the previous section, magnetic reconnection results in shrinking of 
magnetic loops between two converging magnetic fluxes of opposite polarity and 
the charged particles entrained over the magnetic loops are accelerated (see Figure 
1). In this section we discuss the results of numerical studies that confirm these 
predictions. 
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In what follows, we consider different domains of magnetic reconnection which 
were modeled by solving the isothermal MHD equations numerically in a uniform 
mesh using a Godunov-type scheme [ 22] [23] 1241 . 

In order to integrate the test particle trajectories a data cube obtained from the 
MHD models is frozen in time and then 10,000 test particles are injected in the 
domain with random initial positions and directions and with an initial thermal dis- 
tribution. For each particle the relativistic equation of motion is solved 

^(ymu) ^q(E + uxB), (10) 
dt 

where m, q and u are the particle mass, electric charge and velocity, respectively, 
E and B are the electric and magnetic fields, respectively, 7 = (l — u^/c^^ ' is the 
Lorentz factor, and c is the speed of light. The electric field E is taken from the 
MHD simulations 

E= vxB + riJ, (11) 

where v is the plasma velocity, 7 = V x B is the current density, and rj is the Ohmic 
resistivity coefficient. The resistive term above can be neglected because its effect 
on particle acceleration is negligible ||231 . These studies do not include the particle 
energy losses, so that particles gain or loose energy only through the interactions 
with the moving magnetized plasma. 

We note that since we are focusing on the acceleration process only, very sim- 
ple domains can be considered which represent only small periodic boxes of entire 
magnetic reconnection or turbulent sites. For this reason, the typical crossing time 
through the box of an injected thermal particle is very small and it has to re-enter the 
computational domain several times before gaining significant energy by multiple 
scatterings. Thus, whenever a particle reaches the box boundary it re-enters in the 
other side to continue scattering 1.24.1 . 



3.1 Particle Acceleration in 2D domains 

Figure |2]presents an evolved 2D MHD configuration with eight Harris current sheets 
in a periodic box ll23l (see also [ 12 |). The initial density profile is such that the total 
(gas plus magnetic) pressure is uniform. Random weak velocity fluctuations were 
imposed to this environment in order to enable spontaneous reconnection events and 
the development of the magnetic islands. 

Figure|2]clearly shows the merging of islands in some locations and the resulting 
stretching or shrinking which provides appropriate conditions for particle accelera- 
tion. KGLl 1 |23| find that an increase of the parallel velocity component is mostly 
observed within shrinking islands and in current sheets (see the red and yellow zones 
in Figure |2]i, while the increase of the perpendicular component is observed mostly 
near and within stretching islands and between current sheets (see the green and yel- 
low zones in Figure|2|i. This complex behavior is related to the degree of island de- 
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Fig. 2 Topology of the magnetic field represented as a gray texture with semi-transparent color 
maps representing locations where the parallel and perpendicular particle velocity components are 
accelerated for a 2D model with B- = 0.0 at time 6.0 in the code units. The red and green col- 
ors coiTespond to regions where either parallel or perpendicular acceleration occurs, respectively, 
while the yellow color shows locations where both types of acceleration occur The parallel com- 
ponent increases in the contracting islands and in the current sheets as well, while the perpendicu- 
lar component increases mostly in the regions between current sheets. White boxes show regions 
that are more carefully analyzed in |23 | paper The simulation was performed with the resolution 
8192x4096. 10,000 test particles were injected in this snapshot with the initial thermal distribution 
with a temperature corresponding to the sound speed of the MHD model. (From 123J.) 



formation and the particle direction and speed. Within contracting magnetic islands 
or current sheets the particles accelerate predominantly through the first order Fermi 
process, as previously described, while outside of the current sheets and islands the 
particles experience mostly drift acceleration due to magnetic fields gradients [23 j. 
In Figure |3] the first of these effects is zoomed in an example of a single test proton 
which is trapped in a shrinking island and is accelerated. Its parallel speed increases 
while the gyro rotation slows down. This results in an exponential growth of the 
kinetic energy of the particle (as shown in the right panel). 

Similar results were found in 2D collisionless PIC simulations fl^, This 
implies that the first-order Fermi acceleration process within shrinking islands is 
not restricted to collisionless physics or kinetic effects as previously suggested and 
described by PIC simulations (e.g. ifTOl [121 [131 and references therein). This accel- 
eration mechanism in reconnection sites works also in collisional plasmas, under 
the MHD approximation, as shown above and, in fact, MHD codes present an easier 
way to study the physics of particle acceleration numerically. 
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Fig. 3 The case of a shrinking island where the particles accelerate efficiently. It con'esponds 
to region Rl of Fig. |2] The left panel shows the trajectory of a single test proton trapped in this 
contracting island. We see two small magnetic islands on both sides of the central elongated island 
which are merging with it. This process results in the contraction of the central island. The right 
panel shows the exponential increase of the particle energy. The proton orbiting around the center 
of the magnetic island increases its energy increment after each orbit. (Extracted from |23'i|.) 

3.2 Acceleration Near and Within Current Sheets 

In the current sheet zones (regions R4 to R6 of Figure |2) we can also identify a first 
order Fermi acceleration due to simple particle scattering between the converging 
flows entering both sides of the current sheet (or even in merging/shrinking islands 
which are just forming there), as described in yj. 

In zones above and below flie current sheets particles possibly experience pre- 
dominantly a drift acceleration driven by non-uniformities of the magnetic field 
(e g j [4rj). Generally, this effect is less efficient than the first order Fermi process 
in merging/contracting islands and results in smaller acceleration rates. The origin 
of this effect is due to the net work done on a charge by the Lorentz force (Eq.fTTTi in 
a zone of non-uniform large scale magnetic field. The principal equation governing 
this is the scalar product of the particle velocity (or momentum) and the accelera- 
tion by the convective electric field, ~qv x B. In uniform magnetic fields, the energy 
gain and loss acquired during a gyroperiod exactly cancel, so in result no net work 
is done, AW — 0. 

Figure|4]zooms in the details of the acceleration of a test particle near and within 
a single (Sweet-Parker shaped; see also Section 3.4) current sheet. Before the par- 
ticle reaches the current sheet discontinuity it is drifted by the plasma inflow and 
the increasing gradient of B as it approaches the current sheet. When it enters the 
discontinuity (the white part of the trajectory in the left panel), it bounces back and 
forth several times and gains energy (which increases exponentially as shown in 
the right panel of Figure |4]i due to head-on collisions with the converging flow, on 
both sides of the magnetic discontinuity, in a first order Fermi process, as described 
in GL05 degouveiaOS. At the same time it drifts along the magnetic lines which 
eventually allow it to escape from the acceleration region. Therefore, we see two 
mechanisms: a drift acceleration (dominating outside of the current sheet) and first 
order Fermi acceleration inside the current sheet. These processes naturally depend 
on the initial particle gyroradius, since it determines the amount of time the particle 
remains in the acceleration zone before escaping. 
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Fig. 4 The case of acceleration near and within a single cun'ent sheet with a Sweet-Parker config- 
uration like region (R6 in Fig.|2). The left panel shows the trajectory of a test proton approaching 
the diffusion region. The color of the trajectory corresponds to the particle energy (which increases 
from red to yellow and then finally to white when the particle reaches the current sheet). The 
right panel shows the evolution of the particle energy. In the model of Sweet-Parker reconnection 
presented in this figure we used explicit large resistivity coefficient T) = 10^^ in order to make 
reconnection fast. The grid size in the model was set lo Ax = 1/ 1024. (Extracted from |23].) 

3.3 2D versus 3D simulations 

The results presented in the previous sections were obtained for 2D models without 
a guide field. This means that in this case the magnetic lines creating the islands are 
closed and a charged particle can be trapped indefinitely in such an island. The pres- 
ence of a guide field normal to the plane of Figure |2] opens the magnetic loops and 
allows the charged particles to travel freely in the out-of-plane direction. Moreover, 
the islands evolve much slower in the presence of a strong guide field. 

Figure |5] depicts the time evolution of the kinetic energy of the particles which 
have their parallel and perpendicular (red and blue points, respectively) velocity 
components accelerated for three models of reconnection. The kinetic energy is nor- 
malized by the proton rest mass value, i.e., it is actually 7— 1) that is plotted, where 
7 is the Lorentz factor In the 2D model without a guide field (as in the models 
studied in the previous section) there is an exponential growth of energy mostly due 
to the acceleration of the parallel component which stops after the energy reaches 
values of 10^-10^. From that level on, particles accelerate their perpendicular com- 
ponent only with smaller linear rate in a log-log diagram. In the 2D model with a 
weak guide field B-=Q. 1 normal to the plane of Figured there is also an exponential 
acceleration of the parallel velocity component, but due to the presence of the weak 
guide field, this component accelerates further to higher energies at a similar rate as 
the perpendicular one. This implies that the presence of a guide field removes the 
restriction seen in the 2D model without a guide field and allows the particles to in- 
crease their parallel velocity components as they travel along the guide field, in open 
loops rather than in confined 2D islands. This result is reassured by the 3D model 
in Figure |5] where no guide field is necessary as the MHD domain is fully three- 
dimensional. In this case, we clearly see a continuous increase of both components, 
which suggests that the particle acceleration behavior changes significantly when 
3D effects are considered, i.e. where open loops replace the closed 2D reconnecting 
islands. 
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Fig. 5 Kinetic energy evolution of a group of 10* protons in 2D models of reconnection with a 
guide field strength B;=0.0 and 0. 1 (top and middle panels, respectively). In the bottom panel a fully 
3D model with initial B,=0.0 is presented. The colors show how the parallel (red) and perpendicular 
(blue) components of the particle velocities increase with time. The contours correspond to values 
0.1 and 0.6 of the maximum number of particles for the parallel and perpendicular accelerations, 
respectively. The energy is normalized by the rest proton mass energy. The background magnetized 
flow with multiple cuiTent sheet layers is at time 4.0 in Alfven time units (t = L/V^, where L is the 
size of the computational domain and Va is the Alfven speed corresponding to the initial magnetic 
field in the system) for all models. (From |23 1.) 

With the parametrization considered, the gyroradius of a proton becomes com- 
parable to the size of the box domain when its Lorentz factor reaches a value of a 
few times 10*. The largest islands in the system can have sizes of a few tenths of the 
size of the box. These rough estimates help us to understand the energy evolution in 
Figure |5] In the case with no guide field (top panel of Fig.|5]l, the exponential par- 
allel acceleration stops right before the energy value 10* is reached. After this, the 
rate of acceleration significantly decreases. This occurs because the Larmor radius 
of the particles has become larger than the sizes of biggest islands. Therefore, from 
this level on the particles cannot be confined anymore within the islands and the first 
order Fermi acceleration ceases. After that, there is a much slower drift acceleration 
(of the perpendicular component only) caused by the gradients of the large scale 
magnetic fields and acceleration between islands. If a guide field is inserted in such 
a system (as in the model of the middle panel of Fig.|5]l, the picture is very similar. 
However, since the particles are now able to travel along the guide field, their paral- 
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lei velocity component also continues to increase after the 10"^ threshold. Of course, 
in the 3D model, the particles follow the same trend (bottom panel of Fig.|5]l. 

While in two dimensional MHD models without a guide field the parallel accel- 
eration saturates at some level, in the presence of an out-of -plane guide field or in 
three dimensional models this saturation effect is removed. 



3.4 Acceleration in 3D Sweet-Parker Reconnection 

In the Sweet-Parker model of reconnection of two large scale magnetic fluxes of 
opposite polarity |50, 44], the speed of reconnection, i.e. the speed at which two in- 
flowing magnetic field lines annihilate by ohmic dissipation, is roughly rj/A, where 
A is the width of the current sheet discontinuity (Figure 1) and Tj is the Ohmic mag- 
netic resistivity. The entrained plasma follows the local field lines and exits through 
the edges of the current sheet at roughly the Alfven speed, Va=B/ (4;rp)'''^, where 
p is the local density. Thus using momentum flux conservation it is easy to demon- 
strate that the resulting reconnection speed is a tiny fraction of the Alfven speed, or 
Vrec ~ Va5^'/^, where S — LVa/t] is the Lundquist number and L is the length of 
the current sheet. Due to the typically huge astrophysical sizes of the reconnection 
sites, S is also huge for Ohmic diffusivity values (e.g., for the interstellar medium, 
S ~ 10^^) and this makes the Sweet-Parker reconnection very slow. However, ob- 
servations require a reconnection speed close to Va in several circumstances (e.g., 
in solar flares). A way to speed up reconnection is to invoke plasma instabilities, as 
for instance, the stream instability which makes resistivity anomalously large in the 
relation above |45|. Another way is to consider the presence of turbulence in the 
current sheet [28], a process that will be described in 13 .5 1 

In the model shown in the top of Figure |6] KGL12 t24ll investigated the acceler- 
ation of thousands of particles in a Sweet-Parker current sheet but, in order to make 
reconnection fast, they employed a diffusivity coefficient r\ ~ 10^^ expressed in 
code units which, due to the numerical diffusivity, is several orders of magnitude 
larger than the typical Ohmic diffusivity in astrophysical environments and besides, 
makes the Sweet-Parker reconnection in the simulation efficient. The time evolution 
of the energy distribution for the accelerating particles is shown for this model in 
the top left panel of Figure |6] Initially, the perpendicular acceleration dominates, 
because the volume in which particles are injected is much larger than the current 
sheet. The perpendicular acceleration, due to a drift of the magnetic flux, starts be- 
fore the particles reach the reconnection region [23 1. The distribution of particles 
does not change significantly until t = 1.0. Then, a rapid increase in energy by 
roughly four orders of magnitude appears for a fraction of particles. We observe a 
big gap between the energy levels before and after these acceleration events, which 
is also evident in the particle energy spectrum depicted in the subplot of the same 
diagram. The events are spread in time because particles gain substantial energy at 
different instants when crossing the current sheet. The energy growth during this 
stage is exponential. This is clearly due to the first-order Fermi acceleration pro- 
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cess, as stressed before [|4]|23]ll4] and already shown in Figure |4] We note that, as 
in Figure |5] the particles accelerate at smaller rates after reaching the energy level 
~ 10^, because the thickness of the acceleration region becomes smaller than their 
Larmor radii. 

Although the Sweet-Parker model with an artificially enhanced resistivity results 
in a predominantly first-order Fermi acceleration, only a small fraction of the in- 
jected particles is trapped and efficiently accelerated in the current sheet (see the 
energy spectrum of the accelerated particles in the bottom right of Figure |6]l. This is 
because the acceleration zone is very thin. 




Fig. 6 Left column: Particle kinetic energy distributions for 10,000 protons injected in the Sweet- 
Parker reconnection (top), fast magnetic reconnection (middle), and purely turbulent (bottom) do- 
mains. The colors indicate which velocity component is accelerated (red or blue for parallel or 
perpendicular, respectively). The energy is nonnalized by the rest proton energy. Subplots show 
the particle energy distributions at / = 5.0. Right column: XY cuts through the domain at Z = of 
the absolute value of current density | J| overlapped with the magnetic vectors for the Sweet-Parker 
reconnection (top), fast reconnection (middle), and purely turbulent domains (bottom). For the top 
and middle models with large scale current sheets it was employed 6o; = 0.1, 7J = lO^"', and a 
resolution 256x512x256, while for the bottom pure turbulent model it was employed Bq, = 0.2 
and a resolution 128x256x128. (From [241 .) 
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3.5 Acceleration in 3D reconnection sites with Turbulence 



As remarked before, Lazarian & Vishniac ll28l (see also the Chapter by Lazarian in 
this book) proposed a model for fast reconnection that does not depend on the mag- 
netic diffusivity (see also ifTsl ). Given the fact that MHD turbulence is ubiquitous 
in astrophysical environments, this may be a universal trigger of fast reconnection. 
The predictions of this model have been tested successfully by numerical simula- 
tions II22I l24l which confirmed that the reconnection speed is of the order of the 
Alfven speed and independent of resistivity. An important consequence of fast re- 
connection by turbulent magnetic fields is the formation of a thick volume filled 
with small scale magnetic fluctuations. In order to test the acceleration of particles 
within such a domain, KGL12 1 24| introduced turbulence within a current sheet with 
a Sweet-Parker configuration (as described in the previous paragraph) and followed 
the trajectories of 10,000 protons injected in this domain. 

The middle left panel of Figure |6] shows the evolution of the kinetic energy of 
the particles in this case. After injection, a large fraction of test particles accelerates 
and the particle energy growth occurs earlier than in the Sweet-Parker case (see also 
the energy spectrum at f = 5 in the detail at the bottom right of the same diagram). 
This is explained by a combination of two effects: the presence of a large number of 
converging small scale current sheets and the broadening of the acceleration region 
due to the turbulence. Here, we do not observe the gap seen in the Sweet-Parker 
reconnection, because particles are continually accelerated by encounters with sev- 
eral small and intermediate scale current sheets randomly distributed in the thick 
volume. The acceleration process is clearly still a first order Fermi process, as in the 
Sweet-Parker case, but more efficient as it involves larger number of particles, since 
the size of the acceleration zone and the number of scatterers have been naturally 
increased by the presence of turbulence. 



3.6 Acceleration by Reconnection in Pure 3D Turbulent 
Environments 

The bottom left panel of Figure|6]shows the kinetic energy evolution of accelerated 
particles in a domain with turbulence only, i.e., without large scale magnetic flux 
tubes and thus no large scale current sheet. This could be the situation in typical 
diffuse MHD environments like the interstellar, the intracluster and intergalactic 
media. One of the fundamental points of the Lazarian & Vishniac theory |28 | is 
the fact that whenever there is MHD turbulence, there will be fast reconnection of 
the turbulent magnetic field lines from the injection to the dissipation scales of the 
turbulence ll33l l47l l34l . Therefore, particles will be able to accelerate while trapped 
within these multiple current sheets at all scales. 

We see in Figure |6] that the acceleration is less efficient at the beginning and a 
much smaller fraction of particles is accelerated than when a large scale current 
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sheet is present, as in the middle panel of Figure |6] In the later case, the converging 
flow on both sides of the large scale current sheet brings approaching scattering 
centers that undergo only head-on collisions with the particles allowing a continuous 
growth of the particle energy until the saturation level. In pure turbulence, however, 
the absence of a large scale converging flow results in a random particle scattering 
into both approaching and receding small scale magnetic fluctuations (although at a 
smaller rate), so that the overall acceleration is a second-order Fermi process. 



4 Discussion and Conclusions 

Advances both in the understanding of magnetic reconnection in the MHD regime 
and improvement on high energy observations have lately motivated the studies 
of particle acceleration in reconnection sites of astrophysical sources and environ- 
ments. 

In this Chapter, we reviewed particle acceleration in 2D and 3D (collisional) 
MHD domains of magnetic reconnection. It has been shown that particles can be 
efficiently accelerated by reconnection through a first order Fermi process within 
large scale current sheets (build up by large scale converging magnetic fluxes), spe- 
cially when local turbulence is present. The later makes the reconnection fast f28l 
and the volume of the accelerating zone thick Il23l l24l. The particles trapped within 
the current sheet suffer several head-on scatterings with the contracting magnetic 
fluctuations as originally predicted by GL05 [4] (see also 1 1 1) and undergo an expo- 
nential growth in their kinetic energy, as demonstrated numerically by KGLl 1 ||231 
and KGLl 2 |[24l . In a Sweet-Parker configuration where the reconnection speed is 
made artificially large, the acceleration rate is slightly smaller because of the thin- 
ner current sheet, but it is also a first-order Fermi process. In contrast, in pure 3D 
turbulent environments (with no large scale current sheets), particles with gyroradii 
smaller than the injection scale of the turbulence are accelerated through a second 
order Fermi process while interacting with both approaching and receding small 
scale turbulent current sheets. This process can be particularly important for cosmic 
ray acceleration in diffuse turbulent environments like the interstellar, intracluster 
and intergalactic media, while the first order Fermi acceleration in large scale cur- 
rent sheets can be relevant particularly in stellar coronae, compact sources (like the 
accretion disk coronae in AGNs, microquasars, etc.; see Figure 7 as an illustrative 
example), and highly magnetized flows (like AGN, microquasar and GRB jets). 

It has been also shown that the acceleration within reconnection sites (especially 
in the presence of turbulence) works both, in collisional and collisionless environ- 
ments. The acceleration by magnetic reconnection in the 2D MHD regime ll23l l24l 
successfully reproduces the results obtained with more complex (collisionless) 2D 
PIC codes (e.g. |fT2l[T3l). This proved that the acceleration in reconnection regions 
is a universal process which is not determined by details of the plasma physics or 
kinetic effects. However, in the collisional case, only the injected particles with Lar- 
mor radii near the MHD scales are effectively accelerated. This injection problem 
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Fig. 7 From left to right the figure shows: the HST image of M87 AGN; a schematic representation 
of the expected magnetic field stmcture around the accretion disk and the central black hole (as in 
[4 1); a schematic representation of the reconnection zone with the two converging magnetic fluxes 
of opposite polarity as in a Sweet-Parker configuration; and a 3D MHD simulation of magnetic 
reconnection with turbulence injected within the current sheet to make reconnection fast (as in 

ISl). 



can be solved using hybrid codes able to resolve both the kinetic and the MHD 
scales. 

It should be noted that ll43l also investigated particle acceleration in MHD recon- 
nection regimes. However, they concluded that MHD should not be a good approxi- 
mation to describe the process of acceleration by reconnection. This is because their 
3D numerical simulations were performed in a fully resistive MHD regime. There- 
fore, they obtained an efficient particle acceleration due to the high electric field 
induced by the resistivity term only (see Eq. 11) and an absorption of most of the 
available magnetic energy by the electrons in a very small fraction of the charac- 
teristic time of the MHD simulation. This led them to conclude that resistive MHD 
codes are unable to represent the full extent of particle acceleration in 3D reconnec- 
tion. KGLl 1 and KGL12 f23l[24ll on the other hand, explored particle acceleration 
in a nearly ideal MHD regime where only small numerical resistivity was present. 
In this case, the contribution of a resistivity induced electric field is negligible when 
compared to the advection component, namely, the electric field resulting from the 
plasma motion in the magnetized medium, E = —qv x B. 

It has been also demonstrated that the acceleration of energetic particles in 2D 
and 3D reconnection domains shows substantial differences 1231 . This calls for fo- 
cusing on realistic 3D geometries of reconnection. The numerical studies l24l 
have also revealed that apart from the first order Fermi acceleration, additional less 
efficient acceleration processes, like drift acceleration due to non-uniform magnetic 
fields and second order Fermi, also interfere in the process. 

All these numerical studies of particle acceleration have neglected the time evo- 
lution of the MHD environment. This is in general valid since this is much longer 
than the particle time scales. In fact, when considering, for instance, the accelera- 
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tion within large scale current sheets with turbulence, particles are accelerated by 
magnetic fluctuations in the turbulent field and interact resonantly with larger and 
larger structures as their energy increases due to the scatterings. In a steady state 
turbulent environment, as considered here, particles will see on average the same 
sort of fluctuation distribution, so that after several Alfven times, one should expect 
no significant changes in the particle spectrum due to the evolution of the large scale 
MHD environment |9|. Nonetheless, this evolution may be important when consid- 
ering more realistic non-steady environments and when calculating real spectra and 
loss effects (e.g. |35, 20 1). It may be also relevant when considering the (second 
order Fermi) acceleration in pure turbulent environments ( as in Figure |6] bottom 
panel). In this case, electric fields arising from slow modes (betatron acceleration) 
can be relevant to the acceleration process making it twice as larger since the be- 
tatron term dB/dt contributes as much as the electric field term (Eq. 11) in the 
second order process. In forthcoming studies when considering more realistic non- 
steady environments, MHD data cubes varying in time should be used specially in 
pure turbulent studies. 

It has been seen that analytical studies of the first order Fermi process in large 
scale current sheets predict that [4, 14]: (i) the acceleration rate is similar to that 
for shock acceleration; and (ii) the energy power law spectrum of the accelerated 
particles can be even harder [14] than the one predicted for shock acceleration and 
independent on the reconnection velocity |4|. These predictions, although based on 
very simplified assumptions can be, in principle, tested with the numerical simula- 
tions. However, a larger parametric space considering, e.g., different ratios between 
the initial Alfven (or reconnection) speed and the light speed, and different ampli- 
tudes of the injected turbulence (which speeds up reconnection within the large scale 
current sheet) must be still performed in order to assess the sensitivity of both the 
acceleration rate and the particle spectrum to the physical conditions in the recon- 
nection domain. The results from the numerical simulations so far |f 23l l24l indicate 
that the acceleration rate is nearly independent of the initial Alfven (and the recon- 
nection) speed, at least for non-relativistic MHD plasmas |9|. Besides, depending 
on the particle energy, this rate can be larger than the one predicted by shock ac- 
celeration. The determination of the dependence of the energy spectral index and 
the fraction of accelerated particles at high energies with the initial and boundary 
conditions is more complex (see Figure |6]l and requires further numerical studies, 
particularly considering the effects of particle feedback. Recently |13| have eval- 
uated the energy spectrum considering merging islands in 2D collisionless recon- 
nection and obtained a hard energy power-law index, —1.5. However, 3D effects 
and particle feedback have to be accounted as well, and this is a study in progress 
(e.g, |9|). Also, forthcoming studies will require the inclusion of the relevant loss 
mechanisms of the accelerated particles, in order to assess the importance of the 
acceleration by reconnection in comparison to other processes (e.g., shock acceler- 
ation) and to reproduce the observed light curves of the sources. Preliminary studies 
in this direction are also in progress ll20l . 

The numerical simulations shown here focused on proton acceleration. Although 
applicable to electrons too, the numerical integration of the electron trajectories is 
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much longer. Nonetheless, these tests are also needed. In particular, it has been 
also suggested that in an electron-positron pair plasma the pairs could annihilate in 
compressed reconnection sheets [14], so that this could be relevant to acceleration 
by reconnection, e.g., in pulsar winds and relativistic jets in general, as stressed 
before. 




Fig. 8 2D relativistic-MHD jet. Left diagrams: top panel depicts the topology of the magnetic field 
represented as gray texture for a 2D relativistic jet at ? = 80.0 in code units with initial unifonu 
longitudinal magnetic field corresponding to a ratio between the gas pressure and the magnetic 
pressure j3 = 1/300, a density ratio between the jet and the ambient medium 0.01; a jet Mach 
number My = 6.0, and a Lorentz factor F = 10.0. The left middle and bottom panels show the 
same diagram but with superimposed semi-transparent color maps representing locations where 
the parallel and perpendicular particle velocity components are accelerated. The red and green col- 
ors correspond to regions where either parallel or perpendicular acceleration occurs, respectively, 
while the yellow color shows locations where both types of acceleration occur The top left panel 
clearly depicts the large scale features typically seen in 2D jet simulations, i.e., the bow shock at 
the head where the light jet beam impacts supersonically a much denser environment, and internal 
shocks (or knots) all along the beam which are driven by the pinch mode of the Kelvin-Helmholtz 
(K-H) instability. We also identify a more turbulent cocoon enveloping the beam which is formed 
by the mixing of the shocked ambient and jet gas. Magnetic islands can be clearly distinguished in 
this region. The left middle panel depicts only the accelerated particles within the beam and in the 
bow shock region, while the bottom left panel depicts the particles which are accelerated mostly in 
the sun'ounding cocoon. The right diagrams show the acceleration rates: the top one corresponds 
to the middle left panel, i.e., to the acceleration regions within the beam and the bow shock and 
therefore, is dominated by first order Fermi acceleration behind shocks; and the bottom right panel 
corresponds to the bottom left panel, i.e., to the acceleration mostly in magnetic islands and thus is 
dominated by first order Fermi due to magnetic reconnection. The simulation was perfonued with 
a resolution 8192x4096. 10,000 test particles were injected in this snapshot with an initial thermal 
distribution with a temperature corresponding to the sound speed of the relativistic-MHD model. 
(From 1251 ; in preparation). 
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Finally, first order Fermi acceleration by reconnection seems to be powerful 
enough in the solar (or stellar) corona and wind and in the Earth magnetotail where 
this mechanism has been more extensively explored. In order to check whether it is 
also really as powerful in relativistic environments, like the surroundings of black 
holes/accretion disks and jets, as suggested in the literature, further studies of the 
present model in the relativistic MHD (RMHD) regime are also required. As an ex- 
ample. Figure m shows preliminary results of 10,000 test particles injected in a 2D 
relativistic MHD jet system [25]. The colors highlight the regions where particles 
are being accelerated, mostly through first order Fermi, to relativistic energies both 
behind shocks and within magnetic reconnection islands. The comparison between 
the acceleration rate diagrams on the right hand side of this figure clearly indicates 
that both mechanisms are competitive [25 1. These results are encouraging and may 
have rather important consequences on particle acceleration and high energy emis- 
sion processes in microquasars, pulsar winds, AGNs, GRBs and their surroundings 
and, therefore, require more extensive studies. To assess the role of this mechanism 
in modeling flares in the spectrum of compact sources is also in order fli 15] l6l [141 . 
Besides, particles accelerated in domains of reconnection, particularly in pure tur- 
bulent regions, may be available as seed populations for further first order Fermi 
shock or magnetic reconnection acceleration in these different systems. 

To summarize, magnetic reconnection is now recognized as an essential process 
not only in astrophysical sources, but also in turbulent environments which in turn, 
are ubiquitous. In this situation the acceleration of particles by reconnection may 
play a vital role, the importance of which should be evaluated with further research. 
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